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Abstract 

New structure function data at small Bjorken x from HERA are used 

along with next-to-leading order QCD evolution to predict a cross section 

for charged-current interactions of ultrahigh energy neutrinos with nucleons. 

This new result is over twice the size of previous estimates and has important 

implications for cosmic ray experiments now underway as well as for KM3 

arrays (cubic kilometer-scale neutrino telescopes) now in the planning stages. 
PACS numbers: 95.55.Vj, 96.40. Tv, 98.70.Sa, 95.85.Ry 
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Introduction The neutrino-nucleon charged current interaction cross section has a fas- 
cinating energy dependence revealing different physical regimes. At the lowest energies, the 
total cross section goes like the neutrino laboratory energy squared due to the non-relativistic 
final state phase space. In the next regime, the energy dependence becomes linear, as one 
expects from Lorentz covariance, dimensional analysis, and the pointlike 4-Fermi interaction. 
These are sufficient to predict that the cross section goes like Gf 2 s, where s is the invariant 
center of mass energy squared. This increase with energy of the cross section has repeatedly 
been verified and makes the neutrino beam a practical experimental tool. However, the 
linear rise with energy is unsustainable as it would eventually violate s-wave unitarity. The 
4-Fermi model breaks down when s (or the magnitude of the invariant momentum transfer 
Q 2 ) approaches the W-boson mass squared M\y 2 ■ Above that energy, the exchange of a sin- 
gle W-boson on an elementary target predicts a cross section which behaves like cr log( M s 2 ), 
where a = G r* M w 2 m 

This expectation is again revised because the nucleon is a composite, not elementary, 
object made of quarks and gluons. In the ultra high energy limit, s M\y 2 , the total cross 
section is dominated by an integration region of nearly constant energy loss, with momentum 
transfer approximately in the range 0.1M W 2 < Q 2 < M w 2 , and x < 1. In this limit the 
number of quarks with small x thus determines the total cross section. Recent small x 
data from the HI and Zeus collaborations at HERA indicate a rather singular growth in the 
number of quarks carrying small momentum fraction, resulting in a substantial enhancement 
for the physical nucleon cross section compared with the one found for scattering from an 
elementary object. 

In this Letter we will combine new results from the HERA electroproduction experiments 
with the total cross section analysis to predict the ultra-high energy (UHE) neutrino cross 
section at energies of 100 TeV and above. This cross section has important implications 
for cosmic ray experiments now underway and for KM3 arrays in the planning stages, the 
so-called neutrino telescopes. Since cosmic ray photons with energies near 10 15 eV are 
strongly attenuated on the intergalactic photon background, at these energies we are left 
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with neutrinos as the only long lived, light, and electrically neutral elementary particle for 
viewing the UHE universe. Even at energies where photons are still available, neutrino 
astronomy opens a new window on the dense, hidden regions of the cosmos which are 
opaque to photons but virtually transparent to neutrinos [1]. Galactic cores provide good 
examples of such hidden regions. Calculations of AGN sources [2] predict a sizable flux of 
UHE neutrinos in the energy range of 10 14 eV and above. Our predictions for the cross 
section in this range are considerably larger than those of our own previous work (denoted 
MR) [3], and the subsequent work of Reno and Quigg (RQ) [4], which were both based on 
incomplete knowledge of small x physics. The situation has changed recently due to HERA 
electroproduction experiments, where an 800 GeV proton beam and a 30 GeV electron beam 
is used to study the process e + p — > e' + X and extract the quark distribution functions at 
small x. A remarkable rise of the structure functions with decreasing x has been observed. 
This rise strongly influences the UHE neutrino-nucleon cross section. A microscopic strong 
interaction process of fundamentally non-perturbative origin (the quark distributions) thus 
profoundly affects electroweak cosmic rays processes which originate in the most distant 
parts of the universe. 

Small x Evolution and the Neutrino Cross Section The kinematics of one-W exchange 
are standard. The cross section can be written in terms of the scaling variables x = -3-r 
and y = -^7-, and the structure functions , , and , which depend on x and 

Q 2 = sxy. We use the relation 2xT\ = and also ignore terms proportional to JF 3 , the 
difference between the number of quarks and antiquarks, which are negligible in the sea 
dominated region of a; <C 1. Thus, 



G F 2 M W 2 M w 2 /-i /-i (l-y + ±y 2 )FZ N (x,sxy) 

2ir s J^ml X Jouml V (Mwi + xy y 



The denominator of Eq. (1) contains the W propagator ^ 2+J 1 - 2 ^ 2 written in terms of x 
and y. The integration is strongly affected by the singularity which is cut off by M\y 2 ■ 
The lower integration limits give a slight underestimate of the true cross section and simply 
reflect ignorance of the structure function in non-perturbative regions. Our procedure is to 



choose some reasonable Q 2 ow as a cutoff and then test that the results are insensitive to its 
precise value. We have found no significant variation over the range 6.4 < Q 2 ow < 64 GeV 2 . 

To continue we need the parton distributions in order to form = x[u(x) + u(x) + 

d(x) + d(x) + s(x) + s(x) + ...]. Information on these quark distributions can be obtained 
from electroproduction, where F 2 ep = x[^(u(x) + u(x)) + |(o?(x) + d(x)) + + s(x)) + ...]. 

The HI [5] and ZEUS [6] [7] groups at HERA have recently measured F 2 ep over a region of 
10" 3 - 7 < x < 10"- 96 and 10 ' 9 < Q 2 < 10 3 - 7 (GeV 2 ). The important integration region of 
Eq. (1) extends up to Q 2 M w 2 , a region where no direct small x parton data on F 2 ep exist. 
Fortunately, QCD evolution via the well known dGLAP [8] evolution equations, together 
with the available F 2 ep data, enable us to effectively bridge this gap. 

Our procedure is to fit solutions to the next-to-leading order QCD-dGLAP evolution 
equations for the parton distributions to the measured structure function F 2 ep . In the region 
of x < 10 -1 , these equations have a remarkably simple two parameter solution when the 
starting distributions have a power-law form in x. We will present the details elsewhere. 
The result is 

xq(x, Q 2 ) = ^"" p S ^lj , Q2l ex PNn(^)M^)P gg (/i + 1, Q 2 )] , (2) 

where P qg (/i + l, Q 2 ) and P gg (/i + l, Q 2 ) are + 1 moments of the quark-gluon and gluon-gluon 
splitting functions at leading plus next-to-leading order, and the overall normalization A and 
the power, /i, are adjustable parameters. Using Eq. (2), we form the structure function with 
the Ansatz, q(x) = q(x) = u(x) = d(x) = s(x) = 2c(x) = 2b(x) so that P 2 P = ^-xq(x) and 
= 8xq(x). This simple form for the structure function gives a high quality fit to the 
global set of HERA data. Our best fit yields A = 0.011 ± 0.002 and /i = 0.40 ± 0.03; this 
result is displayed together with the HERA data in Fig.l. The X 2 per degree of freedom 
for this fit is 0.85, which is excellent; our quark distributions are in good agreement with 
existing global fits [9] in the regions where such fits have been reported. Our results are 
actually in somewhat better overall agreement with the HERA data than any we have seen 
to date, including those which fit the initial distribution with many parameters. 



With a prediction for Ti{%-, Q 2 ) in hand, we employ Eq. (1) to obtain our cross section. 
In the "large x v region, where the present work has nothing new to add, we use conventional 
multiple parameter fits for the quark distributions, splicing our result with well established 
structure function results [9] for x > 0.2. We find that the cross section for 50 TeV < E v < 
100 TeV is not very sensitive to exactly where the splice is made, and it has virtually no 
effect for E v > 100 TeV. For 50TeV < E v < 50PeV, a compact expression which fits our 
result is 

a ,lN {E v ) = exp[-80.56 + 0.822(ln£ 1/ ) - 0.0231(ln J E iy ) 2 ] , (3) 

with E v in units of TeV. In Fig. 2, we compare this new result with our earlier calculation, 
MR [3], which chose to neglect valence quarks, and also with the calculation of RQ [4], 
which restored this contribution. Valence effects are small for energies above ~100 TeV. 
The region of existing high energy neutrino-nucleon cross section data is indicated, along 
with one very high energy data point recently extracted by the HI group who measured 
a(ep — y v + X) [10]. The lone data point is at approximately 10 times the highest energy 
previously reported and is in good agreement with our new result, thus tying together two 
completely independent aspects of HERA measurements. Varying the position of our x 
splice over the range 0.4 > x > 0.01 changes the low energy portion (< 100 TeV) of our 
result by less than the error bars on the HI data point. Our cross section is roughly a factor 
of 2.2 larger than previous estimates over the range 100 TeV < E v < 10 PeV. 

At higher energies the rate of growth of our cross section is larger because the small x par- 
ton distributions are following a "power law" evolution that self-generates an accelerated Q 2 
evolution. It is important to note that this is an a-priori unknown, non-perturbative bound- 
ary condition for the evolution, entirely consistent with the perturbative QCD formalism 
used to study it. Differences between the new and old physics are entirely in this boundary 
condition. Previous formulas were based on exp(^ln^lnlnQ 2 ) evolution of a parton begin- 
ning at x = 1 and creating a shower (the Green function of dGLAP evolution). We note 
that the question of parton saturation or recombination, which has received much attention 
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at the lower Q 2 regions accessible to the HERA experiments, is a "higher twist" effect (sub- 
dominant to all powers of logarithms) and not relevant at the rather large Q 2 0(Mw 2 ) 
we study. Using our distributions self consistently, the highest energy (smallest x) that can 
safely be considered are E v < 10 3 PeV (x > 10~ 7 ). 

Implications for UHE Neutrino Astronomy High-energy neutrino telescopes are multi- 
purpose detectors that are expected to make important contributions in astrophysics, and 
possibly particle physics. Conventional sources of UHE neutrinos include scattering of ~ 10 2 ° 
eV protons off of 3°K background photons; these events produce charged pions which then 
decay to produce neutrinos. The Fly's-Eye collaboration has reported detailed evidence 
for correlated changes in the spectrum and composition of cosmic rays in the region above 
10 18 eV [11]. Expected point sources within our own galaxy include X-ray binaries such as 
Hercules X-l, Cygnus X-3, and the Crab Nebula. Recently, much attention has been given 
to active galactic nuclei (AGN) as sources for the highest energy cosmic rays and large fluxes 
of UHE neutrinos [12]. Learned and Pakvasa [13] have recently discussed the "new physics" 
signal of tau neutrino oscillations at PeV energies. 

A natural "KM3" scale for UHE neutrino telescopes is 1 km 3 . Several current proto- 
types covering a few percent of this volume are now in various stages of deployment. The 
DUMAND, NESTOR and BAIKAL experiments are located in deep ocean water, deep lake 
water and shallow lake water respectively while the antarctic AMANDA project is situated 
under ice. These experiments detect Cherenkov radiation at optical frequencies. There is 
also a realistic possibility of using coherent radio techniques [14] to enhance detection for 
energies >50-100 TeV. 

As an example of what our UHE cross section implies for KM3 telescopes, we present an 
event rate calculation for 1 km 3 of polar ice based on the latest estimates of UHE neutrino 
flux from AGN sources by Szabo and Protheroe [2]. In general, a larger cross section 
translates into a higher event rate. However, if one is interested in upward going (through 
the earth) neutrinos, then the competing effect of neutrino absorption by the earth must 
be carefully considered. An equivalent solid angle can be calculated which expresses the 
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fraction of the original flux which actually reaches the detector [15]. For our calculation we 
made use of the PREM earth structure model [16] from which an integrated nucleon density 
over possible incident neutrino directions can be obtained. Our result appears as part (a) of 
Fig. 3, where the effect of our larger cross section can be seen by comparison with results 
obtained using the RQ cross section. The energy at which 50% of the upward flux fails to 
reach the detector is about 700 TeV in the case of RQ compared with roughly 180 TeV using 
our results. 

In part (b) of Fig. 3 we show a calculation of the differential event rate, the number 
of electromagnetic cascades per year per TeV per cubic kilometer of ice, induced by muon 
neutrinos. Rates for upward and downward going events are indicated separately and dashed 
versus solid lines indicate results using the two cross sections as in part (a). Due to the 
effect of earth shadowing, the upward rates are always smaller than downward rates with the 
difference becoming more pronounced as cascade energy increases. Interestingly, the upward 
rates for our larger cross section are higher than for the MR or RQ cross sections despite 
the lower flux at the detector due to shadowing. In the competition between attenuated 
flux versus increased interaction probability for the surviving particles, the net effect is an 
increase in rate. 

Rate predictions for a specific experimental situation depend on an effective detection 
volume for the apparatus as a function of the event (cascade) energy, which include array 
geometry, receiver characteristics, signal attenuation in the medium, noise characteristics 
and detection threshholds. Given the bulk event rates shown in Fig. 3, the prospects for a 
KM3-scale array of almost any reasonable design are very encouraging as is the corresponding 
potential for scientific discovery. 
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FIGURES 

FIG. 1. Our best fit (lines) compared with the HI and ZEUS data. Numbers just above each 
line indicate Q 2 values for the data in GeV 2 . For clarity, the data and associated fits are separated 
by integer powers of 3 as Q 2 increases. Thus, we show 3° times .^(^j 8.5), 3 1 times J~2(x, 12— 15), 
and so on thru 3 10 times J-2(x, 5000). 

FIG. 2. Total charged current neutrino-nucleon cross section versus incident neutrino energy. 
The extrapolated 4- Fermi result is shown along with the single high energy HERA data point. Also 
shown is the present ultrahigh energy prediction (solid line) along with previous results of MR and 
RQ. 

FIG. 3. (a) Equivalent solid angle (units of 2ir) for the upward going neutrino flux showing the 
increasing opacity of the earth with energy. Attenuation is greater with the new larger cross section 
(denoted FMR) compared with the RQ cross section estimate, (b) The differential event rate per 
year per TeV versus muon energy for a cubic kilometer of polar ice due to AGN neutrinos. We 
compare results obtained using the FMR and RQ cross sections with separate curves for upward 
versus downward going events. Note that the integrated rate between 10 2 TeV and 10 3 TeV for 
upward going events using the new cross section is 8029 events/year/km 3 . 
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